Electric Field Effect Tuning of Electron-Phonon Coupling in Graphene 
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Gate- modulated low-temperature Raman spectra reveal that the electric field effect (EFE), per- 
vasive in contemporary electronics, has marked impacts on long wavelength optical phonons of 
graphene. The EFE in this two dimensional honeycomb lattice of carbon atoms creates large density 
modulations of carriers with linear dispersion (known as Dirac fermions). Our EFE Raman spectra 
display the interactions of lattice vibrations with these unusual carriers. The changes of phonon 
frequency and line- width demonstrate optically the particle-hole symmetry about the charge-neutral 
Dirac-point. The linear dependence of the phonon frequency on the EFE-modulated Fermi energy 
is explained as the electron-phonon coupling of mass-less Dirac fermions. 
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The interaction between electrons and quantized lat- 
tice vibrations in a solid is one of the most fundamental 
realms of study in condensed matter physics. In par- 
ticular, the electron-phonon interaction in graphene and 
its derivatives plays an important role in understanding 
anomalies of photoemission spectra observed in graphite 
[l| and graphene Q , the non-linear high energy electron 
transport in carbon nanotubes (3, R IE IE , as weu 
as phonon structures in graphite [8j, l9( and carbon nan- 
otubes [1 LIE El- 

Traditionally, electron-phonon interactions are inves- 
tigated through chemical doping, in which the charge 
carrier density is varied by introduction of impurities. 
The electric field effect (EFE) is an alternative method 
for changing the charge carrier density effectively in low- 
dimensional systems. The EFE has proven very success- 
ful in graphene, a single atomic sheet of graphite, where 
unconventional integer quantum Hall effect 

[[111! h as 

revealed physics linked to the uniqueness of the electronic 
band structure near the charge neutral Dirac points (Fig. 
IHa)). 

We measured Raman spectra of optical phonons in 
graphene where large densities of free electrons or free 
holes are modulated by the EFE. We discovered that 
the even parity long wavelength optical phonon (the 
graphene G band) has marked dependence on gate volt- 
age and the induced charge density. The dependence of 
phonon frequency and line-width on the EFE induced 
charge density demonstrates that the intriguing physics 
of mass-less Dirac fermions with particle-hole symmetry 
is encoded in the electron-phonon interaction. 



Raman studies of graphite [14| are at the forefront of 
research on carbon based materials. The recent avail- 
ability of few-layer and single-layer graphene [IE Ll6j ]. 
has stimulated great interest in Raman scattering in 
such novel and exciting systems. For example, di- 
mensional crossover was observed in Raman spectra of 
thin gra phi tic films as a function of multilayer thickness 
[ttL 0, HH . In the work reported here, Raman spec- 
troscopy emerges as an insightful method to probe the 



EFE in a single atomic layer and the phonon dynamics 
that are associated with the two dimensional (2D) Dirac 
fermions. 

We focus on the doubly degenerate optical phonon of 
E<2g symmetry at ^1580 cm -1 , known as the G band. 
We also report on the smaller impact of the EFE on 
the second-order band at ^2700 cm" 1 , known as the D* 
band. These two bands are prominent Raman features 
in graphene 17 , IE Il| • G band is the optical phonon at 
long wavelengths and D* band is associated with a two- 
phonon state in which each phonon has a large wavevec- 
tor. Our experiments demonstrate that the G band is 
markedly sensitive to coupling with Dirac fermion exci- 
tations at small wavevectors (long wavelengths), while 
the reduced impact of the EFE on the D* band offers 
insights on the coupling to particle-hole pairs at large 
wavevectors. 

The interaction of the G phonon with small wavevec- 
tor particle-hole pairs has a crucial dependence on the 
onset energy e c for vertical transitions (zero momentum 
transfer) of an electron from tt valence band to a 7r* con- 
duction band state as shown in Figs. LUb) and [He). The 
7r — > 7r* transitions need to satisfy the requirements of 
the Pauli principle. For this reason, e c is directly linked 
to the position of the Fermi surface as e c = 2\Ep\, where 
\Ep\ is the absolute value of the Fermi energy. The re- 
sponse of the G phonon to the EFE is thus inextricably 
related to the electronic band structure in the vicinity of 
the Dirac point. 

Graphene samples were prepared using the mechanical 
method described in Ref. [2fJ . The sample was mounted 
on degenerately doped p-type silicon with a 300-nm thick 
Si02 as gate dielectric. A gold electrode was thermally 
evaporated on one end of the graphene layer to contact 
it to the electrical ground (Fig. Did)). A gate bias volt- 
age V g was applied to the silicon substrate to induce 
charge carriers in the sample by the EFE. The sample was 
mounted in a variable temperature cryostat with optical 
access where it was cooled by a continuous flow of cold 
Helium gas. The Raman measurements were performed 
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FIG. 1: (color online). Electronic band structure of graphene 
(a-c) and experimental set-up (d). (a) n and n* bands of 
graphene in momentum space. K and K' (Dirac points) 
are the corner points of the first Brillouin zone. For neu- 
tral graphene, the ir band is completely filled shown as blue 
shading in this diagram. Near Dirac points, the band struc- 
ture exhibits isotropic linear dispersion shown as two cones 
at K and K'. (b) shows the region near Dirac points in n-type 
graphene. e c is the onset energy for vertical electron-hole pair 
transitions, (c) is like (b) for p-type graphene. (d) Upper left 
is an optical image of the sample. 



in a backscattering configuration (Fig. [IJd)) with the 
488-nm line of an argon-ion laser as exciting radiation. 
The laser spot size on the sample was ~ 25 fim in diam- 
eter and the power was kept below 5 mW. The scattered 
light was collected into a Spex-1404 0.85-m double grat- 
ing spectrometer and spectra were recorded by a liquid 
nitrogen cooled CCD camera. The spectral resolution in 
these measurements is ^2 cm . 

Typical spectra showing G and D* bands at cryostat 
cold finger temperatures ~10 K are displayed in Fig. [5] 
Remarkably, both peak position and line- width are found 
to be modulated by the gate voltage. The changes in- 
duced by the EFE are much larger for the G mode than 
for the second-order D* band. We also show for compar- 
ison a G band for thick (>10 nm) samples (Fig. Hfb)) 
taken under similar experimental conditions with no gate 
voltage applied, which serves as a bulk sample spectrum. 

The EFE induced charge in the single layer of graphene 
can be modulated by gate voltage according to the rela- 
tion n = C g (V g - V Dirac )/e, where C g , e, and V Dira c 
are the gate capacitance, electron charge and the gate 
voltage corresponding to the charge-neutral Dirac-point. 
The changes in graphene Raman spectra as a function 
of V g demonstrate that the EFE modulated charge den- 
sity results in significant changes in the frequency and 
lifetime of phonons. 

Figure displays the G-mode frequency log and width 
Tq as a function of V g . It is remarkable that uaiVg) and 
Tq(V 9 ) both exhibit nearly symmetric changes relative to 
the value of V g ~ 18 V. This symmetry is naturally linked 
to the symmetry of the electronic band structure that 
occurs at the Dirac point shown in Fig. [2(a), from which 
we determine Voirac — 18±2 V. In particular, u>g(V 3 ) has 



3. 

'</> 
c 




(b) 



graphite G band 
•\ 1-1 OK 

V 



1560 1570 1580 1590 1600 1610 



monolayer graphene 


(c) 


D* band 


T=12K 






V 

SI 


Jr l V 


80V 




40V 








20V 




OV 




% 
-40V 



2660 2680 2700 2720 2740 



Raman shift (cm' ) 

FIG. 2: (color online). Low temperature Raman spectra and 
their Lorentzian fits (smooth overlapping curves), (a) Evolu- 
tion of the G band of a graphene monolayer with gate voltage 
V g . The upper left inset is a schematic representation of car- 
bon atom motion in the G band. The dashed line is a guide 
to the eye. (b) Spectrum of the G band of a thick (>10 nm) 
graphite layer, (c) Evolution of the D* two-phonon band of 
a graphene monolayer with V g . 



a minimum near the Dirac point, revealing that increases 
in charge density of either sign result in stiffening of the G 
mode. On the other hand, Tq sharply decreases as \V g — 
VDirac\ increases, showing that longer phonon lifetimes 
are linked to higher particle/hole density. These results 
offer key insights into the coupling of Dirac fermions to 
lattice vibration modes, as discussed below. 

We first consider Tg shown in Fig. [4ja). Its marked 
reduction (from 15 to 6.5 cm -1 ) at higher carrier densi- 
ties can be understood as a Landau damping of phonons 
in which the mode decays into particle- hole pairs. Ex- 
amples of such processes are shown in the inset to Fig. 
HJa) and in Fig. 0|b). Phonon decay processes are real 
transitions that conserve energy and momentum. The 
small wavevector G phonon can only decay into particle- 
hole pairs represented by vertical transitions that have 
vanishingly small wavevector transfer (Fig. HJb)). When 
charge carriers are induced by EFE in ideal graphene 
at low temperatures, the Fermi level Ef changes as 
Ef{ji) = —sgn(n)hvF\/ {^\n\) , where vf = 10 6 m/sec 
[l2l [13} is the Fermi velocity and n is charge density. 
To satisfy the Pauli principle, Landau damping is al- 
lowed only when \E F \ < ftw G /2 (Figs, gjb) and|4jc)). 
This is the relatively low density regime in which G band 
phonons are damped via decay into vertical particle-hole 
pairs. 

The total change in Tg is written as ATq = 
^ g(V Dirac)— where ^g ^ s ^ ne res idual line- width from 
processes that are not related to Landau damping. The 
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FIG. 3: (color online). G band energy (squares) and G band 
width (circles) extracted from Fig. 2(a). The vertical dotted 
line is the approximate position of the charge neutral Dirac 
point that is estimated from the symmetry of the data. The 
upper scale is obtained with n = C g (V g — VDirac)/e, where 
C 9 =115 aF/^tm 2 [3 EH- 



Landau damping phonon decay rate ATq is the proba- 
bility of creating energetically allowed electron-hole pairs 
via the electron-phonon coupling in unit time. Therefore, 
the electron-phonon coupling strength can be estimated 
from Arc. A calculation based on the Fermi golden rule 
yields 
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where A uc is the area of the graphene unit cell, M is 
carbon atom mass, and D is the electron-phonon cou- 
pling strength. The measured 8.5 cm" 1 Ar^ yields 
D = 14.1 eV/A, which seems to be consistent with the 
value predicted by density functional theory [2~0 ]. 

The results in Fig. [4|a) exhibit that the measured Tq 
has departures from the step-like function (blue dashed 
line) expected for ideal graphene with uniform charge 
density distribution. In graphene samples, however, 
there are presumably local density variations due to self- 
doping [221 ] . chemical adsorbants [15} and charged impu- 
rities trapped in the substrate [23j. Different locations 
of the graphene sheet have thus different local contri- 
butions to damping of G phonons, smoothing out the 
sudden changes at ±fta>G/2. Considering a non-uniform 
spatial variation of ~ ± 3 x 10 11 cm -2 in charge density, 
we obtained a reasonable fit for Tg shown as the red solid 
curve in Fig. Ufa). 

The EFE modulation of wg , that has a minimum near 
the Dirac point, is shown in Fig. Hfd) as a function of 
Fermi energy. For carrier densities above ~ ± 5 x 10 11 
cm~ 2 where non-uniformity is no longer dominant, the 
G-band phonon energy increases linearly with changes in 
Fermi energy. This intriguing result directly links the 




-300 -200 -100 100 200 300 

Fermi energy (meV) 

FIG. 4: (color online). Graphene G band damping (a-c) and 
energy renormalization (d-f). In (a) and (d), blue dashed 
lines and red solid lines are the fits for ideal and non-uniform 
graphene respectively. In all cases, the same electron-phonon 
coupling strength (D = 12.6 eV/A) is assumed. The vertical 
dotted lines are the determinations of the Dirac point shown 
in Fig. The insets are Feynman diagrams for electron- 
phonon coupling applicable to the case of the G phonon. (b) 
represents the broadening of the G phonon due to decay into 
particle-hole pairs, (c) indicates that the G phonon decay into 
electron-hole pair is forbidden by the Pauli principle at high 
EFE induced charge densities, (e) is for the renormalization 
of the G phonon energy by creation of virtual electron-hole 
pairs, (f) shows that virtual electron-hole pair transitions 
with energy ranging from to 2|_Bf| are forbidden. Only the 
diagrams for n-type graphene are shown in panels (b), (c), 
(e), and (f). Diagrams for the p-type cases are similar. 



EFE-modulated loq to 2D Dirac fermion properties, as 
discussed below. 

The coupling between G phonon lattice vibrations and 
Dirac fermions is allowed by graphene lattice symmetry. 
The carriers residing in the honeycomb lattice respond to 
the dynamical perturbation of G mode lattice vibration 
by creation and annihilation of virtual long-wavelength 
electron- hole pairs across the gapless Dirac point, which 
in turn, renormalizes the G phonon energy. The energy 
range of the virtual electron-hole pairs allowed by Pauli 
principle is decided by the position of Fermi level (Figs. 
EJe) and[UJf)). When the graphene is charge neutral, the 
onset energy is zero. If graphene is doped with electrons 
or holes, the onset energy is twice the Fermi energy. The 
difference of the G band energy in neutral and charged 
graphene is thus given by the renormalization of electron- 
hole pairs with energy ranging from to 2\Ep\. 

Carrying out the calculation explicitly with time de- 
pendent perturbation theory 2J] and the linear disper- 
sion near the Dirac point, we find that the change of 
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G-band phonon energy is described by 



^ G ln| 2 !f F !~^ G |} (2) 



2\E F \ + ftw G 



where u> G is lug at Dirac point, A = 



m. In 



2-rrhuJGMv'j, 

the high density limit, the linear term dominates over the 
logarithm term. A fit of the linear segments in Fig. 0Jd) 
(blue dashed lines) with the linear term in Eq. (2) gives 
D = 12.6 eV/A, close to the value 14.1 eV/A determined 
above from the analysis of the ojq data. The red smooth 
line in Fig. IDJd) is the fit that takes into account the 
impact of non-uniformity described above. 

It is worth noting that = ^f 3 - «154 meV, which 
depends only on one parameter. Our experimental result 
gives 195 meV. 

Finally, we consider the impact of the EFE on the 
D* Raman mode shown in Fig. Etc). Although both 
peak position and width of the G and D* bands exhibit 
similar gate dependence behavior, the magnitudes of the 
changes seen in Fig. Etc) are only ~10% of those of the G 
band. This smaller gate response can be explained by the 
impact of electron-phonon interactions for the two large 
wavevector phonons that contribute to the second-order 
Raman D* band. The renormalization of the phonon fre- 
quency due to electron-phonon interactions is inversely 
proportional to the energy of particle-hole transitions at 
the phonon wavevector [24J. In this context, the weaker 
response of the D* band to the EFE suggests that rele- 
vant EFE-modulated particle-hole transitions here have 
energies about an order of magnitude larger than those 
in the renormalization of the G band [26j . Further quan- 
titative description could be carried out by the analysis 
of the doubly-resonant Raman processes in the D* band 

0- 

In conclusion, we observed EFE-tunable electron- 
phonon coupling in Raman spectra of single layer 
graphene. The gate voltage dependences of phonon fre- 
quency and damping reveal charge-tunable interactions 
of optical phonons with Dirac fcrmion transitions across 
a vanishing bandgap. These results uncover physics that 
links the EFE with electron-phonon coupling and demon- 
strate venues to probe fundamental interactions and to 
characterize atomic scale structures and devices. 

Additional Remark: During the preparation of this 
manuscript, we became aware of related work on simi- 



lar systems from other groups [27], |28|, [29J, |30| • 
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